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At any resolution level of wavelet expansions the physical observable of the kinetic energy is 
represented by an infinite matrix which is "canonically" chosen as the projection of the operator 
—A/2 onto the subspace of the given resolution. It is shown, that this canonical choice is not 
optimal, as the regular grid of the basis set introduces an artificial consequence of periodicity, and it 
is only a particular member of possible operator representations. We present an explicit method of 
preparing a near optimal kinetic energy matrix which leads to more appropriate results in numerical 
wavelet based calculations. This construction works even in those cases, where the usual definition 
is unusable (i.e., the derivative of the basis functions does not exist). It is also shown, that building 
an effective kinetic energy matrix is equivalent to the renormalization of the kinetic energy by a 
momentum dependent effective mass compensating for artificial periodicity effects. 



I. INTRODUCTION 

Wavelets are commonly used for analyzing and for a 
compact storage of complex distributions like two di- 
mensional images, temporal signals, and even for solving 
partial differential equations. Goedecker and Ivanov [l| 
solved the Poisson equation, Cho et al. employed wavelets 
in solving the Schrodinger equation for Hydrogen-like 
atoms [2]. Using this tool, all electron calculations were 
also performed within the framework of the local den- 
sity approximation applying pseudopotentials and su- 
percells Car-Parinello algorithm [3| and in Ref. [5] 
a new approach for magnetic ordering was presented. 
Arias and his coworkers developed Kohn-Sham equations 
based wavelet method |6j] , and also tested for various sys- 
tems (e.g., and [1]). 

Wavelet analysis is a popular label applied for the con- 
cept of multiresolution analysis (MRA), which covers a 
systematically refined basis function set of Hilbert spaces. 
We refer to basic textbooks (see e.g., (91,1101) for the de- 
tails. 

In our previous works we have shown that the sur- 
roundings of a molecule can be described at a rather 
rough resolution level [ill We have also demon- 

strated [l^l that electron-electron cusp singularity of the 
two-electron density operator can be easily reproduced by 
the method of multiresolution analysis. In [14| we have 
studied the detail structure of the wave function at var- 
ious refinement levels using MRA. An adaptive method 
was also developed for identifying the fine structure lo- 
calization regions, where further refinement of the wave 
function is necessary without solving the eigenvalue equa- 



tion in the whole subspace expanded by the basis func- 
tions of the given resolution level. 



II. SYSTEMATIC ERROR IN FINITE 
RESOLUTION EIGENFUNCTIONS 

While having studied the question, which physical re- 
gions of the potential need a high resolution expansion of 
the wave function, we have solved numerically the matrix 
form of the one particle Schrodinger equation of exactly 
solvable models, with the Hamiltonian 



H 



1 



V{x). 



(1) 



The algebraic representation of the Schrodinger equation 



(2) 



for the ith excited state is derived by considering that 
according to the MRA construction, at the resolution 
level M , the Hilbert space TL is approximated by one of 
its subspaces 7i'*^l = spaTi{sMi{x)\l G Z}, where the 
orthonormal basis functions 



SMi{x) = 2^-^/2^(2^'^ a; - £) 



(3) 



are the scaled and translated versions of the "mother" 
scaling function s{x). The translated scaling functions 
are "sitting" on an equidistant grid of grid length 2^*^. 
The series of subspaces H'*^' {M co) "approximates" 
in a given sense the complete Hilbert space Ti, and the 



2 



projectors 



P. 



M 



\sMt){sMl\ 



(4) 



of "approximate" the identity operator (although, 

in the strict mathematical sense the limit of Pm does not 
exist). 

By inserting the approximate identity Pm into ([2]) one 
arrives at 



(5) 



Multiplying ([5|) by the basis element sjv/j from the left 
results in 



[sMj\H\sMl) {sMl\^i) = Ei{sMj\^i 



(6) 



How well approximation ([5]) works is far from being un- 
derstood. Nevertheless, the above method of algebraiza- 
tion is conventional, and later on, we will refer to this 
procedure as "canonical". 

Of course, the canonical method includes the solution 
of the eigenvalue problem 



[M] 



CM£ 



(7) 



of the Hamiltonian matrix — {sMj\H\sMi) ■ The 

eigenvalue -E]*^' is only an approximation to the exact 
eigenvalue Ei (an upper bound for the ground state), 
and the eigenvectors cmi define an approximation 



lei 



(8) 



of the wave function '^i{x). One can not expect, of 
course, that ([8]) gives a better result than the best ap- 
proximation 

,[M] 



SMt\^i) SMI 



(9) 



in the subspace T^I^^l. 

For an illustration, we have chosen the simplest analyt- 
ically solvable model of the potential box. The alterna- 
tive of the free electron problem was singled out, as the 
wave functions should be square integrable in order to 
be able to successfully describe it with matrix methods. 
Fig. [T] shows the exact excited state 4*5, its projection 



5°' to the subspace of resolution M — 0, and the so- 



lution $5°' related to the eigenvalue problem ([7]). The 
Hamiltonian matrix was calculated using the compactly 
supported 6 parameter Daubechies scaling functions 9]. 
As the first derivative of these basis functions exists, the 
kinetic energy matrix elements were determined by 



[M] 



^kA«) = ^ ((-iV)sMj|(-iV)sM£) 



{x)s'Mt{x)dx. 



(10) 
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FIG. 1: Exact (solid line) and approximate (dashed line) 
wave functions 'I's and t&g''' of the potential box with 
L — 15 a.u., W = 100 a.u. The difference of the exact wave 
function and the solution of the canonical eigenvalue problem 
is also shown. For a reference the difference of the exact wave 
function and its projection to subspace Ti'"' is plotted as well. 
Atomic units were used. 



Here we have used the fact, that the momentum operator 
— zV is self adjoint. The potential energy matrix elements 
were calculated numerically with the potential function 



V{x)^ 



if \x\ < L, 
W if > L. 



(11) 



By a careless analysis of the results one can easily draw 
erroneous conclusions. One might think, that at the re- 
gions, where the difference of the exact and approximate 
solution is large, a further refinement of the basis set 
is necessary. Theoretically, this could be accomplished 
by adding wavelets sitting in the regions of large errors. 
Wavelets are localized basis functions of the orthogonal 
complement subspace Wl^l of Ti'^^l in the embedding 
subspace Wl^^+^l = ^[^1 ® Wl*^!. 

According to Fig. [1] the large error regions are located 
at the steepest parts of the oscillating wave function. At 
these places, however, the scaling function expansion can 
not be of bad quality, considering that any linear function 
of the form ax + b can be exactly expanded in T^I*^! at 
any resolution level M . As the wave functions at the 
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steepest parts are almost linear, we do not expect large 
errors in the scaling function expansion. This prediction 
is justified in the third plot of Fig.[Tl where the difference 
of the exact and projected wave function is shown. The 
large deviations in the approximate wave function 
should have a different origin. 

A careful study of the first part of Fig. [T] leads to the 
conclusion, that the solution of the eigenvalue problem 
has (apart from small irregularities) an oscillatory form 
similar to the exact solution. The essential difference be- 
tween the two wave functions is that the approximate 
wave function has a slightly smaller wavelength 

than the exact one, leading, of course, to a larger ki- 
netic energy. In the case of the excited state i = 5, e.g., 

(^'5|r|*5) = 0.1352, while (^ifllTl^if') = 0.1389. The 
same general experience was gained by studying other 
excited states of the box model as well as those of the 
harmonic oscillator. 

In the following sections we will analyze the reasons, 
why the above effects appear, and suggest a possible so- 
lution. 



III. THE KINETIC ENERGY OPERATOR IN 
THE SUBSPACES n^"'^ 



in a different Hilbert space. However, to avoid any weird 
representations, we would like to keep some properties of 
the "canonical" approach, which are summarized below. 

1. Considering, that the scaling functions are usually 
(but not necessarily) real, according to pil)) 



[M] 



T, 



[M] 



(12) 



2. Using the definition of sm£{x) and after a simple 
variable transformation in pop . one arrives at the 
shift invariance property 



T. 



[M] 



[M] 



(13) 



indicating that the kinetic energy is represented by 
a band matrix. 

3. If the scaling functions are compactly supported on 
the interval [0, N — 1] (as in the case of Daubechics- 
N bases), 

t\ ^ — - J s'j^,j(>{x)s\,jf^{x)dx — 0, 



if 1^1 > N -2. 



(14) 



As in the previous example the systematic error in the 
wave function occurred in regions without a potential en- 
ergy contribution, we conclude, that the effect is due to 
the representation of the kinetic energy operator. Pos- 
ing this question is not as heretical as one might think 
for the first sight. In the original formulation of matrix 
mechanics by Heisenberg, Born and Jordan, the physical 
quantities are represented by infinite matrices satisfying 
the appropriate commutation rules. There is no specific 
prescription for the determination of their matrix ele- 
ments. On the other hand, Schrodinger works in the 
Hilbert space of the square integrable functions Ti., with 
a specific prescription for the operator representation of 
physical quantities, in particular, T = —A/2. Von Neu- 
mann has shown [1^ the equivalence of both descriptions 
with the abstract Hilbert space equipped with linear op- 
erators for physical quantities. One can consider here 
that the subspace T^I^^l C 7i is itself an infinite dimen- 
sional separable Hilbert space, and as such, can serve for 
a complete description of any quantum mechanical sys- 
tem, on its own right. The significant difference of the 
multiresolution expansion from the usual atomic orbital 
expansions is that the latter span a finite dimensional 
subspace, which is in principle unable to describe a quan- 
tum system in all details. 

As the subspace Ti^*^! essentially differs from the com- 
plete space of square integrable functions, it is natural 
that the elements of the infinite matrix if correspond- 
ing to the physical quantity of the kinetic energy differs 
from the matrix elements rj*^' = {sMj|— A/2 jsMf), since 
—A/2 is the operator representation of the kinetic energy 



Here iV is the number of parameters which define 
the mother scaling function. 

4. Using definition ^ in PH)) leads to a scaling prop- 
erty 



rp[M] _ „2Mrp[a] 



(15) 



5. In three spatial dimensions a direct product basis 
function set \j1j2j3) = SAfji (a;i)sMj2 (a;2)sMj3 (^3) 
is used. As the Laplacian is a simple sum of second 
derivatives according to the three spatial variables, 
orthonormality of the basis functions results in 

(J1J2J3I- A/2 14^2^3) = TlfilSn^.S,,,, 

+Sj,i,6j^i^T^jffl. (16) 

Consequently, the case of the 3D kinetic energy op- 
erator is straightforwardly reduced to the one di- 
mensional representation. 



Hermiticity (|12p and translational invariance (jl3p are 
natural requirements for any operator representations of 
the kinetic energy. Though we are not obliged to keep 
property (HU, this is one of the most attractive features 
of using compactly supported basis sets, thus this pre- 
scription is applied as well. Transformation of the ki- 
netic energy matrix elements with increasing resolution 
(like the scaling property (fT5|) ) will be discussed later. 
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For the matrix i^T^*'^! representing the kinetic energy we 
require the foUowings at any resolution level M 



K 



[M] 



A' 



[M] 



^[Af] ^ ^[Af] ^ ^[Af] 



[M] _ 



0, if 1^1 > TV -2. 



(17) 
(18) 
(19) 



In three dimensions we additionally apply 



K 



[M] 



(20) 



The role of operators assigned to observables is to pro- 
vide the possible and expectation values of the corre- 
sponding physical quantities. A proper representation of 
the kinetic energy should give the known fc^/2 eigenvalues 
with eigenvectors which give reasonable approximations 
of the free electron wave function e**^^. In the following 
considerations we will prove that the best approximation 
Pmg^'^^ is really an eigenvector of any matrix satisfy- 
ing the requirements p7 P "(|19 p . The question remains to 
clarify how well the equality 



i^ji {sMi\e ) - 



[SMj 



1 



(21) 



is satisfied. 

Defining the Fourier transform by /(^) — 

using definition ^ we 

have 

{sMAe''^) - 2-*^/2e'fc-^(2^)i/2 s(-fcM), (22) 

where kM = 2~^^fc is the scaled wave number. Using this 
expression in the left hand side of ([2T|l 



E7^[Af], I ikx\ 



(23) 



It is clear that just the shift invariance (|T8|) ensures that 
the projection of the free electron wave function to T^I*^! 
is an eigenfunction of ifl^^l, with the eigenvalue 



(24) 



eel, 

As according to PT|) if]*^' is Hermitian, its eigenvalues 
£[*^l(fc) are real and 

AM] tkMi 

■I ^ 



i.e., e[*^l(fc) is symmetric in k. The above natural phys- 
ical requirements are satisfied for all reasonably chosen 

The fundamental question is how well e[^^l(fc) approx- 
imates the free electron kinetic energy fc^/2. In order to 
understand this, some properties of e^^'^\k) will be stud- 
ied below. As the argument k^.j in definition ex- 
ponentially decreases with increasing resolution M, it is 
natural to consider the Taylor expansion of el^l (fc). Due 
to the symmetry of e[*^l(fc) all its odd order derivatives 
should be zero in /c = 0. This condition is equivalent to 
Efez ^^"+^A:f = for any n = 0, 1, . . ., which follows 
from ITSl) . In the ideal case e[^^l(fc) would be equal to 
fc^/2 requiring 



Am 



(0) 



[M] 



0, (26) 



,-2Af ^ ^2 T^m 

eez 

^^(0) = (-l)"2-2"^^5]£2"xfl 



(0) = -2-2A^^^2^fl ^ 1, (27) 
" " 0, (28) 



for n > 2. In Appendix El we have shown that for the 
canonical kinetic energy matrix sum rules (jAl[) and (|A3p 
ensure that and (P7)) are satisfied. It is easy to see, 
however, that regardless of the choice of Kf^^ the equiva- 
lence of (k) and k'^/2 can never be perfect. According 
to (HH) £M(/c) is a finite trigonometric polynomial, con- 
sequently it is 2tt periodic, and can not coincide with the 
free electron energy in the whole range ~oo < k < oo. 

Fig.[2]shows that the canonical kinetic energy function 
approaches the "Brillouin zone" boundaries at A; = —it 
and fc = TT with a horizontal tangent. This difficulty can 
never be resolved by choosing any set of k'^/'^^ values for 
2n periodic e^*^! (k) functions, which leads to the conclu- 
sion, that any multiresolution quantum mechanics loses 
its applicability for the energies with |fc| = 2tt/X « tt. 
In other words, for the wavelength of the particle we get 
the condition A/2 <^ 1. It is not surprising from the 
physical point of view, since 1 is the grid length of the 
scaling function basis set at resolution level M = and 
no wave functions with a wavelength comparable to the 
grid length are expected to be described in a satisfactory 
manner. If the resolution increases, however, the scaling 
property (jlSp results in 



E 

^ez 



rl^^e-^'^-'' = 22*^e[l(2-*'^fc). (29) 



(25) 



Since the argument of eian decreases exponentially with 
increasing resolution, the quality of the Taylor expan- 
sion becomes increasingly better, as it can be traced in 
Fig. [21 in the close neighborhood oi k = 0. Although it is 
very satisfying that the canonical kinetic energy can re- 
produce the exact values in the infinite resolution limit, 
in a practical calculation, however, one should stay at 
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A. The Fourier series approach 
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FIG. 2: The canonical kinetic energy function ecan(fc) (dashed 
Une) of the Daubechies-6 basis set compared to the free parti- 
cle kinetic energy A;^/2 (solid line). The kinetic energy e^'^ik) 
of the Fourier type approximation with the same number of 
matrix elements is plotted by a dotted line. Atomic units 
were used. 



a relatively low resolution M. At these resolutions the 
function eitS (fc) performs rather poorly for larger fc, and 

there is a reasonable hope to find matrix elements Kf^^ 
which provide significantly better approximation of the 
function 



IV. ASSIGNMENT STRATEGIES FOR THE 
OPTIMAL KINETIC ENERGY MATRIX 

There are several possible decisions for choosing the 
matrix elements k}-^^^ in order to approximate by a 
finite trigonometrical polynomial expansion of the form 
As a first remark, we recall the scaling property ([15]) 
of the canonical kinetic energy matrix. Since the values 
are not determined by a mathematical formula 
(similar to PU)) ). their scaling behavior can not be de- 
rived. On the other hand, however, the necessary require- 
ment (|27l) implies that in the leading order the kinetic 

energy matrix elements should scale as K^"^^ ~ 2^*^. We 
would like to emphasize, that choosing the scaling for- 
mula 



(30) 



is by no means a must, we have decided to apply it in 
order to decrease the many degrees of freedom of the 
problem. After this, we still have to determine the values 
of the zeroth level matrix elements, for which we consider 
two different philosophies. 



The best possibility we can expect using an expression 
of the form ((24)) is that the function fc^/2 is correctly 
described in the interval (— tt, tt). This, however, requires 
the infinite Fourier series expansion 



2 " 6 +^ 



cos{£k). 



(31) 



Applying (|T9l) leads to a truncation of ((3T|) at £ = iV — 
2. It is clear that identifying Kf^ with the expansion 
coefficients of the truncated series would not satisfy any 
of the criteria (|26 p - (l^51) . It is an elementary requirement 
that a particle with zero momentum should have zero 
kinetic energy (criterion ([26]) ). On the other hand, in the 
large resolution limit we should recover the exact kinetic 
energy, and as we have discussed above, this is equivalent 
to ([27]) . Considering these arguments, we suggest the 
following truncation process 



■N 



1 (ill 

2 \^ 6 



for ^ = 0, 

for 1 < 1^1 < iV - 3, 
for \l\= N - 2. 

(32) 

automatically satisfies 
sum rule (|26p , whereas with an appropriate choice of the 
normalization factor a, ([27]) can also be fulfilled. The ki- 



The definition kI"'^ = 2"^^^^, 



■N 



netic energy function (fc) calculated according to ((24|) 
using the matrix elements determined by (j32p is plotted 
in Fig. [2l According to the figure, the Fourier type ap- 
proach results in a weaker quality approximation than 
the canonical calculation, especially in the low energy 
region. This effect can be traced in Tab. [H both the to- 
tal energy and the wave function deviate more from the 
exact quantities than those of the canonical calculation, 
except for the wave function of higher excited states. 



B. The Taylor series approach 

As we have seen, the quality of the £M(fc) for smaller 
k values is essential both in rough resolutions and also in 
the limit M ^ oo. This leads to the conclusion that the 
Taylor expansion of the kinetic energy function should 
satisfy as much conditions of (PS)) - (|^5|) as possible with 
the given number of non-zero matrix elements i^T^*^' . We 
suggest the following scheme for determining the optimal 
kinetic energy matrix. According to (fTH]) and (fTOj) the 
number of non-zero, essentially different matrix elements 
is iV — 1, offering too much freedom in the optimization 
process. Consequently, we have decided to keep only one 
independent parameter t and to define the kinetic energy 
matrix elements by 



K 



[M] 



(33) 
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TABLE I: The error of the total energy -B]"''^ and the wave function "1>[''''^ of the Fourier method compared to the ones of 
the canonical quantities Ef'^ and for the Daubechies-6 basis set. Ei and are the exact total energy and wave function, 
respectively. Label i = 1 indicates the ground state, whereas i = 2, 3, 4, 5 are the successive excited states of the potential box 
(fTTj). with L = 15 a.u., W = 100 a.u. 
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with 



1 for £ = 0, 
t for 1^1 = 1, 
te{t) for 2 < 1^1 < iV - 



(34) 



Quantities ti{t) are determined by the solution of the 
hnear system of equations 



/ t2 \ 



( 1 


1 


1 


\ 


24 


34 . 


[N - 2)4 




26 


36 . 


{N - 2f 




22(^-3) 


32(JV-3) 


. (7V-2)2(^- 


-3) J 



-1/2-t \ 



(35) 



The normalization factor 



Af-2 



-2i-2^^2^f(i). 



(36) 



It is easy to verify that the conditions (p6|) and pT)) are 
always satisfied by these values and ([28|) fulfills until n = 
TV -3. 

Notice, that equations ([55)1 have a solution even in the 
case of iV = 4, i.e., for the Daubechies-4 scaling functions, 
where the canonical kinetic energy matrix elements are 
not defined at all, as the application of formula (fTO|) re- 
quires the derivative of the scaling function, which does 
not exist in this case. For an illustration we have carried 
out a calculation for the potential 



if -L < a: < 0, 
V{x) = {Vq if < x < L, 
W if \x\ > L, 



(37) 



with Daubechies-4 scaling functions and with the Taylor 
series based method outlined above. The optimal value of 



the parameter t was determined as described later. Fig. [3] 
shows the exact and the M = level approximated wave 
function of the 2nd excited state of the model (|37p. This 
state was selected in order to demonstrate the capabili- 
ties of the method, not only in the free electron case but 
in the classically unavailable regions (right hand side of 
the box, Ei < Vq), as well. It is seen that the construc- 
tion of the kinetic energy matrix by the Taylor expansion 
method leads to a rather satisfactory result even at this 
low resolution approximation. 

For finding the optimal values of t in the construction 
([34l) - (l36|) we have plotted the deviation of the approx- 
imate total energy from the exact one as a function of 
the parameter t in Fig. |4l Similarly, the norm of the dif- 
ference of the exact and approximated wave function of 




FIG. 3: Exact (solid line) and approximate (dashed line) wave 
functions 'I'l and of the potential model (|37p . with 

L = 15 a.u., W — 100 a.u. and Vq = 0.5 a.u., using the 
Daubechies-4 scaling function in the Taylor series method. 
Atomic units were applied. 
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FIG. 4: The deviation of the approximate total energy E[ 
and wave function <l>f'''^ from their exact counterparts Ei and 
$1 (solid lines) of the potential model (|37p . with L = 17 a.u., 
W — 10 a.u. and Vb = 0.01 a.u., using the Taylor series 
method. As references, horizontal lines indicate the values 
of the canonical (dashed line) and the Fourier series method 
(dotted line) results. Sign Q stands for the results of the 
Daubechies-4 scaling function, x , □ and A mean Daubechies- 
6, Daubechies-8 and Daubechies-10 calculations, respectively. 
Atomic units were applied. 



TABLE 11: The recommended kinetic energy matrix parame- 
ter t for the Taylor series based construction (|33p - p6p in case 
of Daubechies-A^ scaling function sets {N — 4, 6, 8, 10). 



N 
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10 
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-0.54 


-0.47 


-0.57 


-0.58 



parameters L, W and Vq of model ((37)) we could realize, 
that the optimum t values are essentially potential in- 
dependent. Extending the scope of the investigations to 
finer resolution levels Af = 0, 1, 2, 3, 4 we have found that 
the position of the minima of the error curves stabilizes at 
M — I, there is no significant change at fine resolutions. 
Further, we have studied the exactly solvable model of 
the harmonic oscillator with the potential 

V{x) = ^x^. (38) 

We have experienced the same behavior as in the case 
of the previous physical system (|57|) . The positions of 
the optimum t values are independent of the excitation 
level i, of the resolution M and also of the value of the 
potential width uj^^. Finally, we conclude, that to a 
good approximation, the optimum values of parameter t 
in the kinetic energy matrix definition can be 

chosen model independently, which works also for many 
low lying excited states. In Tab. [TT] we summarize our 
recommendations for the the best values of t for various 
Daubechies basis function sets. As the error curves are 
rather flat around the minima, the values in the table 
can be changed by ±10 percent without a considerable 
change in the quality of the total energy and wave func- 
tion approximations. 



the model potential ([57]) is also shown. It can be realized, 
that in quite a broad range of the parameter i, the kinetic 
energy matrix of the Taylor's expansion method leads 
to a considerably better quality result than the canoni- 
cal choice, regarding both the total energy and the wave 
function. The Fourier series method, however, leads to 
unacceptable solutions. 

According to Fig. |4l a minimum in the error functions 
can be identified which defines the optimum value of pa- 
rameter t. We would like to emphasize, that around the 
minima the curves are fiat, even for the energy deviations 
(consider, the logarithmic scale on the vertical axis). The 
error curves are similar to all of the lower energy excited 
states, and even the positions of the minima are almost 
the same. As the worst example, we adduce the case of 
Daubechies 8 scaling function where the average value of 
the optimal t calculated for the energy and wave func- 
tion error curves of several excited states is 0.57, with a 
standard deviation of ±0.06. Also, by playing around the 



V. SUMMARY 

Studying exactly solvable models in the framework of 
multiresolution analysis we have found that the repre- 
sentation of the kinetic energy operator plays an essen- 
tial role in the quality of the results achieved by approxi- 
mate solutions at a given resolution level M. The regular 
grid of the scaling function basis set introduces an arti- 
ficial consequence of periodicity. Instead of a free par- 
ticle, the MRA expansion describes rather an electron 
with a momentum dependent effective mass m*{k) = 
fc^/(2e[^-^l (fc)), where the function e[^'^l(fc) is determined 
by the matrix elements of the kinetic energy matrix. We 
have shown that in the case of resolution level M — 
the MRA expansion loses its applicability if the kinetic 
energy approaches or exceeds the value £^kin = 7r^/2 a.u. 
(corresponding to the limit value k = ir). By increasing 
the resolution, the applicability range extends exponen- 
tially as fc < 2^ IT, Ekin < 2^^^7r^/2 a.u., and artificial 
periodicity effects disappear in this limit. 
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However, in the numerical practice the level of resolu- 
tion should be kept as low as possible, in order to avoid 
the need for extensive computational resources. We have 
demonstrated that for low resolutions the kinetic energy 
of the numerical calculations is overestimated compared 
to the exact values. The effect is due to the fact that the 
free particle energy is improperly reproduced by its 

27r-periodic approximation el^nik). This reproduction 
can considerably be improved by introducing alternative 
matrix elements of the kinetic energy matrix instead of 
the canonical ones. Both the total energy and wave func- 
tion improvements are well pronounced at low resolutions 
as it can be traced from Fig. |4l A close optimal, system 
and eigenstate independent choice of the kinetic energy 
matrix elements is derived from the formulas of the Tay- 
lor series approach and from our suggestion for 
its parameter value t in Tab. [TTl At an arbitrary res- 
olution level M the kinetic energy matrix elements are 
calculated as 



condition is satisfied for the Daubechies basis sets with 6 
or more parameters. 

Proof of HAl]) . At any resolution level M the scaling func- 
tion basis set is capable to exactly expand any constant 
function , consequently, for any x 



^4 = 1 



(A4) 



with the expansion coefficients c[*^' = 2 Differenti- 
ating, multiplying by s'j^jq{x)/2 and integrating one gets 



Proof of ItA^} . The basis set {sMe\i G exactly expands 
the identity function, i.e., for all x 



K 



[M] 



,[M] 



E4 *A«(a;) = X 



(A6) 



xsMiix)dx ^2-^^^/^{fii+£), (A7) 



where TV determines the number of essential matrix ele- with the appropriate expansion coefficients 
ments, as 77^ „, = if I j - £1 > TV - 2. Tab. |TTT] lists the 

values of T/^ for various Daubechies- scaling function 
sets. In three spatial dimensions expression pO|) should 
be applied. With the suggested method it is possible to 
define a kinetic energy matrix even in the case of the 
Daubechies-4 basis set, where the scaling function is not 
differentiable, thus the canonical approach is not appli- 
cable. 
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APPENDIX A: SOME ELEMENTARY 
PROPERTIES OF THE CANONICAL KINETIC 
ENERGY MATRIX 

We will prove here, that the canonical kinetic energy 
matrix elements defined by (|10p satisfy simple sum rules 
as 



where we have applied definition ([3]), a proper integral 
variable transformation, and the fact that J s{y)dy = 1 
9]. The quantity /ii — J ys{y)dy is the first momentum 
of the mother scaling function. Differentiating (IA6p . mul- 
tiplying by s'j^jq{x)/2 and integrating we arrive at 



^ f SMo{x)dx = i [sA/o(a;)]!°, 



(A8) 



As the scaling functions are square integrable, 
smo{±co) — 0, and 

[A/]^[M] _ 



^ E4"'T^. 



eel 



2-3M/2^^ J2 t]*^! + 2-3*^/2 J2 i TP . (A9) 



Err = 0, 
E^rf" = 0, 



~,2M 



(Al) 
(A2) 
(A3) 



Considering (jAip immediately follows (|A2[) . 
Proof of liA3\) . For the compactly supported scaling func- 
tions of Daubechies with 6 or more parameters the func- 
tion x'^ is still among the exactly expandable functions. 



Em I \ 2 



(AlO) 



with the notation introduced in (1131) . Of course, these 
relations hold only if the canonical kinetic energy matrix 
exists, i.e., if the scaling function is differentiable. This 



where the expansion coefficients are 
4*^1= f x''sMeix)dx ^2-^^'/\ti2 + 2£tii+f), (All) 
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TABLE III: The recommended kinetic energy matrix elements for the Taylor series based construction in case of Daubechies-A'' 
scaling function sets. 





= 4 


TV = 6 


A = 8 


TV = 10 


-'o 


1.3157894737 


1.0269360269 


1.4668325041 


1.5177613012 


n-N 


—0.7105263158 


—0.4826599327 


—0.8360945274 


—0.8803015547 


-'2 


0.0526315789 


-0.0586700337 


0.1207733653 


0.1495444216 


T-JV 
''3 




0.0326358826 


-0.0206082682 


-0.0344316374 


T-JV 
-'4 




0.0047739298 


0.0027102582 


0.0074722170 


T-JV 
-'5 






-0.0002005664 


-0.0013201227 


T-JV 

ie, 






0.0000034864 


0.0001689233 










-0.0000133612 


T-JV 
-'8 








0.0000004634 



after similar steps applied in the previous proof. The no- 
tation J y^s(jj)dy = ^2 was introduced. Differentiating 
(PIOl) leads to 

- j SMo{x)dx^ -2-^^/"^. (A12) 



The second equality follows from partial integration, and 
from the fact that the integrated part disappears due to 
SAfo(±oo) = 0. Substituting (|Alip into (|A12p and using 
and ([X2|) gives (IA3I). 
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